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Semi-realistic nucleon-nucleon interactions applicable to the self-consistent mean-field (both 
Hartree-Fock and Hartree-Fock-Bogolyubov) calculations are developed, by modifying the M3Y 
interaction. The modification is made so as to reproduce binding energies and rms matter radii of 
doubly magic nuclei, single-particle levels in ^"^Pb, and even-odd mass differences of the Sn isotopes. 
We find parameter-sets with and without the tensor force. The new interactions are further checked 
QQ ' by the saturation properties of the uniform nuclear matter including the Landau-Migdal parameters. 

, By the mean-field calculations, interaction-dependence of the neutron drip line is investigated for 

' the O, Ca and Ni isotopes, and of the single-particle energies for the A'^ = 16, 32, 50, 82 and Z — 50 

fS) ' nuclei. Results of the semi-realistic interactions including the tensor force are in fair agreement with 

available experimental data for all of these properties. 

o 

Q ■ PACS numbers: 21.30.Fe, 21.60.Jz, 21.10.Dr, 21.10.Pc 



I. INTRODUCTION 



Mean-field (MF) theories provide us with a good first approximation to the nuclear structure problems. They are 
able to describe the saturation and the shell structure simultaneously, both of which are basic nuclear properties, 
based on effective nucleon-nucleon (NN) interactions. As far as we constrain to the non-relativistic approaches, most 
of the MF calculations have been performed with the Skyrme interaction [H . Finite-range interactions have rarely 
\ been applied, except the Gogny interaction which has the Gaussian form for the central channels. Most popular 
parameter-sets of the Skyrme and the Gogny interactions have been adjusted mainly to the nuclear properties around 
\^ ' the /?-stability. However, it is a question whether such phenomenological effective interactions work well for nuclei 
ly-j far off the ^stability. For instance, whereas role of the tensor force in nucleus-dependence of the MF has attracted 
P\| interest the tensor force is usually ignored in those parameter-sets. 

] Although direct application of the bare NN (and NNN) interaction to the nuclear structure problems [1, H, 0| 
is yet limited to light nuclei or made to medium-mass nuclei but with limited accuracy, guide from microscopic 
00 theories will be valuable even in heavy-mass nuclei. The Michigan 3-range Yukawa (M3Y) interaction which 
was derived by fitting the Yukawa functions to Brueckner's G-matrix, has been used in nuclear structure as well 
^ as in low-energy nuclear reaction studies. There have been a few attempts applying the M3Y-type interaction to 
• ^ MF calculations ^ld\ . In Ref . [lOj , the author has developed an M3Y-type interaction which is applicable to the 
Hartree-Fock (HF) calculations. The original M3Y interaction is incapable of reproducing the saturation and the 
spin-orbit splitting within the MF regime. To cure this problem a density-dependent contact term has been added 
and some of the strength parameters have been modified. Such interactions, which have originally been derived 
from microscopic theories but arc slightly modified from phenomenological standpoints, may be called semi-realistic 
interactions. It has been shown [lO|, that semi-realistic NN interactions could give different shell structure from 
the widely used Skyrme and Gogny interactions. However, the pairing properties have not been taken into account 
in the parameter-set M3Y-P2 that was proposed in Ref. [l^. This implies that the singlet-even channel in M3Y-P2 
is not quite appropriate as long as the pairing interaction is taken to be consistent with the HF interaction, while 
this problem seems to be masked in the HF approximation. In order to apply semi-realistic interactions to the MF 
studies extensively, we explore new parameter-sets of the M3Y-type interaction, taking the pairing pro per ties into 
account. Special attention is paid also to role of the tensor force. Using the recently developed algorithm ITsl. IT^ . 
we apply the new semi-realistic interactions to the Hartree-Fock-Bogolyubov (HFB) as well as to the HF calculations 
of medium- to heavy-mass spherical nuclei. 
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II. M3Y-TYPE INTERACTION 



We express a non-relativistic nuclear effective Hamiltonian by 

2M 



H^=K + Vn; K = Y,^, VN=Y.v^,, (f) 

i<j 



with i and j representing the indices of individual nuclcons. For the effective NN interaction Vij, we consider the 
following form, 



V, 



(C) 



(LS) 



,,(TN) 



(DD) 
V ■ 



n 

(f E)PSE • W.)f'^' + t(™'PTE • [P(r,)]"'™') '^(r,) , (2) 



where = - r^, = |ry |, p.y = (p,j - Pj)/2, Ly = x p^^-, S'.y = 4 [3(s., • ■ fy) - ■ Sj] with fy = ry /ri^-, 

and p(r) denotes the nucleon density. The Yukawa function /n(r) — e^'^"''//i„r is assumed for all channels except 
^(DD) ^YiQ M3Y-type interactions. The projection operators on the singlet-even (SE), triplet-even (TE), singlet-odd 
(SO) and triplet-odd (TO) two-particle states are defined as 

^ 1 - f + ^ 1 + P^ 1 - P. 

2 2 ' 2 2 ' 

_ 1 - p. 1 - p. _ i + p. f+p. 

Pso-^ (3) 

where Pq- (Pt) expresses the spin (isospin) exchange operator. 

We shall start from the M3Y-Paris interaction [3| , which will be denoted by M3Y-P0 in this article as in Ref. ■ 
We change none of the range parameters /x„ of M3Y-P0 in v^'^\ ^(lS) g^j^j^ In M3Y-P0, the longest range 

part in w'^'-'^ is kept identical to the central channels of the one-pion exchange potential (OPEP), Wopep- ^^^^ 
maintain this reasonable assumption. As is well known, the spin-orbit (is) splitting plays a significant role in the 
nuclear shell structure. Even though higher-order effects may account for the observed £s splitting fw], it is desired 
to enhance v^-^^^ in order to describe the shell structure within the MF regime. We here use an overall enhancement 
factor to v^-^^\ which is determined from the single-particle (s.p.) spectrum of ^°*Pb, as will be shown in Sec. IIVI 
Influence of the tensor force on s.p. energies is a current topic, which could be relevant to the new magic numbers in 
unstable nuclei [3i]. We develop two parameter-sets having ■y*-"^^-' without any modification from M3Y-P0, as well as 
a parameter-set in which we impose v^™'' — 0. 

The saturation properties are important to describe many nuclei in a wide mass range. Since it is still hard to 
describe accurately the saturation properties by the bare NN (and NNN) interaction despite certain progress 0], 
it will be appropriate to modify realistic effective interaction so as to reproduce the saturation properties. Density- 
dependence in the effective interaction has been known to be essential in obtaining the saturation. We therefore 
add a density-dependent contact force [l^. The parameter a^"""^-* in v^^^\ power to p, is taken to be 1/3, by 

which the incompressibility JC becomes close to a reasonable value as shown later. On the other hand, Qf(^^) is not 
quite sensitive to K., because the major source of the saturation lies in the TE channel, not in the SE channel [13] ■ 
Although we simply assumed ~ 1/3 in M3Y-P2 11 Oil, this assumption makes it difficult to reproduce 

pairing properties and to avoid instability of the neutron matter |10l. IT^ simultaneously. To overcome this problem, 
we adopt a*^^^^ = 1 in the new parameter-sets. The difference between a*^^^-* and a^™^ may be attributed to the 
difference in origin of the p-dependence; the short-range repulsion in the bare NN interaction in the SE channel while 
primarily the tensor force in the TE channel. Validity of the choice a^^^^ = 1 is further discussed in Sec. IIIII 

The remaining parameters are t„ in w*^^) (except those of w|j^^gp) and tp in v^^^\ We fit them to the measured 
binding energies of ^^O and ^°^Pb, in the HF approximation (see Sec. lIVp . The proton and neutron Fermi energies of 
^"^^Pb, which are primarily relevant to the symmetry energy, are checked additionally. To determine il^^'' and ip^^\ 
we also use the even-odd mass differences of the Sn isotopes, by comparing results of the Hartree-Fock-Bogolyubov 
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TABLE L Parameters of M3Y-type interactions. 



parameters 




M3Y-P0 


M3Y-P2 


M3Y-P3 


M3Y-P4 


M3Y-P5 




(fm) 


0.25 


0.25 


0.25 


0.25 


0.25 


,(SE) 
'•1 


(MeV) 


11466. 


8027. 


8027. 


8027. 


8027. 


(TE) 
'■1 


(MeV) 


13967. 


6080. 


7130. 


5503. 


5576. 


.(SO) 
'•1 


(MeV) 


-1418. 


-11900. 


-1418. 


-12000. 


-1418. 


.(TO) 
''1 


(MeV) 


11345. 


3800. 


11345. 


3700. 


11345. 




(fm) 


0.40 


0.40 


0.40 


0.40 


0.40 


,(SE) 
'■2 


(MeV) 


-3556. 


-2880. 


-2637. 


-2637. 


-2650. 


(TE) 


(MeV) 


-4594. 


-4266. 


-4594. 


-4183. 


-4170. 


,(SO) 
'-2 


(MeV) 


950. 


2730. 


950. 


4500. 


2880. 


.(TO) 
'-2 


(MeV) 


-1900. 


-780. 


-1900. 


-1000. 


-1780. 




(fm) 


1.414 


1.414 


1.414 


1.414 


1.414 


,(SE) 

'■s 


(MeV) 


-10.463 


-10.463 


-10.463 


-10.463 


-10.463 


(TE) 
^3 


(MeV) 


-10.463 


-10.463 


-10.463 


-10.463 


-10.463 


,(SO) 

'•s 


(MeV) 


31.389 


31.389 


31.389 


31.389 


31.389 


,(TO) 
''3 


(MeV) 


3.488 


3.488 


3.488 


3.488 


3.488 




(fm) 


0.25 


0.25 


0.25 


0.25 


0.25 


(LSE) 
'•1 


(MeV) 


-5101. 


-9181.8 


-10712.1 


-8671.7 


-11222.2 


,(LSO) 
'-1 


(MeV) 


-1897. 


-3414.6 


-3983.7 


-3224.9 


-4173.4 




(fm) 


0.40 


0.40 


0.40 


0.40 


0.40 


(LSE) 
'•2 


(MeV) 


-337. 


-606.6 


-707.7 


-572.9 


-741.4 


(LSD) 
''2 


(MeV) 


-632. 


-1137.6 


-1327.2 


-1074.4 


-1390.4 




(fm) 


0.40 


0.40 


0.40 


0.40 


0.40 


,(TNE) 
'■I 


(MeV-fm-2) 


-1096. 


-131.52 


-1096. 


0. 


-1096. 


,(TNO) 
'•1 


(MeV-fm-2) 


244. 


29.28 


244. 


0. 


244. 




(fm) 


0.70 


0.70 


0.70 


0.70 


0.70 


,(TNE) 
'■2 


(MeV-fm-2) 


-30.9 


-3.708 


-30.9 


0. 


-30.9 


,(TNO) 
^•2 


(MeV-fm-2) 


15.6 


1.872 


15.6 


0. 


15.6 


q(SE) 






1/3 


1 


1 


1 


,(SE) 
[■p 


(MeV-fm^) 


0. 


181. 


220. 


248. 


126. 


a(™) 






1/3 


1/3 


1/3 


1/3 


(TE) 


(MeV-fm) 


0. 


1139. 


1198. 


1142. 


1147. 



(MeV-fm) 



(HFB) calculations with the experimental values (see SecIV]). The new parameter-sets of the semi- realistic M3Y-type 
interaction, M3Y-P3 to P5, are tabulated in Table HI For comparison, M3Y-P0 and P2 are also shown. In the set 
M3Y-P3, we keep both wC^^) and the odd-channel (SO and TO) strengths in t;^^') of M3Y-P0. The set M3Y-P4 is 
obtained by assuming v''™'> = 0, while changing tn^^"* and t^'^'^ [n = 1, 2) substantially. In the set M3Y-P5, we 

somewhat modify ^2^'^'' and ^2"^'^' while keeping v^™\ so as to reproduce the binding energies of several doubly magic 
nuclei better than M3Y-P3, as will be shown in Sec. IIVI Thus the number of adjusted parameters are 7, 11 and 9 
(including the overall enhancement factor to u^'"^)) for M3Y-P3, P4 and P5, respectively. It will be useful to compare 
results of these parameter-sets for pinning down which part of the interaction is important to individual physical 
quantities. In particular, role of the tensor force will be of interest. It is remarked that, while schematic tensor forces 
have been introduced into some of the recent MF studies ^, ilSj] , the present M3Y-P3 and P5 is much more 

realistic. 
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III. PROPERTIES OF NUCLEAR MATTER AT AND AROUND SATURATION POINT 



We first view properties of the infinite nuclear matter that are predicted by the semi-realistic NN interactions. In 
the HF approximation, energy of the nuclear matter can be expressed by the following variables: 



P = 



Vs = 



^Pra- , 

err 

^(^Pto- 

ar 

p 

J2^Pra 



E 



CTTpr 



P 



Pp] ^ Ppl + PnT ^ Pnl 



Pp] + Ppl ~ Pn] - Pni 



Ppl ~ Ppl ~ Pn^ + Pnl 
P 



(4) 



p^g. {t = p,n and a =ti J-i which are sometimes substituted by ±1 without confusion) stands for densities depending 
on the spin and the isospin, and is related to the Fermi momentum k-p-ra by 



(5) 



In the interaction of Eq. only w^^^ + t)(^^) contributes to the energy of the uniform nuclear matter. Formulas to 
calculate the nuclear matter energy and its derivatives for given fepro- have been derived in Ref. [l^l- Note that, even 
when superfluidity makes the nuclear matter energy somewhat lower, it is not much different from the energy in the 
HF approximation. 

The spin-saturated symmetric matter is characterized by rjs ~ rjt — "qst — 0, for which we denote feprtr simply by 
/cp- The minimum of the energy per nucleon £ = E/A, given by 



d£ 
'dp 



d£ 



0, 



(6) 



defines the saturation density po (equivalently, /cpo) and energy Sq. The expression |o indicates evaluation at the 
saturation point. As well as po and fo, second derivatives of £ carry basic information of the effective NN interaction. 
Two of the curvatures of £ at the saturation point are called incompressibility and volume symmetry energy. 



jC ^ 1,2 <±±_ 



_ id^£ 



and are related to the Landau-Migdal (LM) parameters /o and /q as 



/C= :^^(l + /o), at 



(7) 



(8) 



where Mq represents the effective mass (fc-mass) at the saturation point. See Ref. [13] for definition of the LM 
parameters. The other curvatures of £ with respect to rjs and ry.st, denoted by and Ost, are defined analogously to 
at and are expressed in terms of the LM parameters go and g'g jlO|. The fc-mass is defined by a derivative of the s.p. 
energy e(k(TT): 



deCkar) 



dk 



kpo 
MS ' 



and is connected to the LM parameter /i by 



M, 





M 



In addition, density-dependence of the symmetry energy, which is represented by a third derivative of £ as 



Ct = ifcp^^ 
2 dkp drit 



d^£ 



t;P 



2 dp drit 



(9) 



(10) 



(11) 
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TABLE II: Nuclear matter properties at the saturation point. 







DIS 


M3Y-P2 


M3Y-P3 


M3Y-P4 


M3Y-P5 


fcpo 


(fm) 


1.342 


1.340 


1.340 


1.340 


1.340 


fo 


(MeV) 


-16.01 


-16.14 


-16.51 


-16.13 


-16.12 


/C 


(MeV) 


202.9 


220.4 


245.8 


235.3 


235.6 






0.697 


0.652 


0.658 


0.665 


0.629 


at 


(MeV) 


31.12 


30.61 


29.75 


28.71 


29.59 


as 


(MeV) 


26.18 


21.19 


20.17 


15.61 


19.56 


ast 


(MeV) 


29.13 


38.19 


36.45 


39.89 


41.01 


Ct 


(MeV) 


22.44 


27.98 


25.30 


17.87 


24.63 



is under interest in relevance to structure of the neutron star crust [19|. These quantities calculated from the new 
semi- realistic interactions are tabulated in Tables HH We here set M = {Mp + Mn)/2, where Mp (Mn) is the measured 
mass of a proton (a neutron) [2^. For comparison, the values obtained by the DIS parameter-set [2l| of the Gogny 
interaction and from M3Y-P2 are also displayed. 

Related to the global systematics of the binding energies and the radii, fcpo ~ 1-33 — 1.34 fm""'^ and Sq ~ —16 MeV 
have been established empirically. Although the M3Y-P3 interaction yields deeper £q than the other interactions, it 
is still within the range of ambiguity in extracting the volume energy from the experimental data (2^ . In practice, 
M3Y-P3 does not yield overbinding for any of the doubly magic nuclei presented in Sec. IIVI 

For the incompressibility, K. « 240 MeV is extracted from the experimental data [2^ . The fc-mass is empirically 
known to be Mq « (0.6 — 0.7)M (23 |. The volume symmetry energy at is important in reproducing global trend of 
the binding energies for the Z ^ N nuclei, and from empirical viewpoints at « 30 MeV seems appropriate [25'|. These 
are fulfilled reasonably well in all the new parameter-sets M3Y-P3 to P5. The choice a^^^-' — 1 contributes to the 
shghtly higher /C in M3Y-P3 to P5 than in DIS and M3Y-P2. 

Global characters of the spin and isospin responses are customarily discussed in terms of the LM parameters. By 
using the formulas given in Ref. [To| . we evaluate the LM parameters for the new semi-realistic interactions, as shown 
in Table Hill It has been known that is small while is relatively large (~ 1) [1^. Although the LM parameters 
should eventually be checked by corresponding excitation modes in actual nuclei, which is beyond the scope of this 
paper, all the semi-realistic M3Y-type interactions seem to have reasonable characters on the spin and isospin channels. 
Not necessarily true for phenomenological interactions such as DIS, this may be linked to the microscopic origin of 
the interactions. In particular, fQ^^EP carries about half oi in the results of the M3Y-type interactions [lo| . 

Figure [1] illustrates S[p) for the spin-saturated symmetric nuclear matter obtained from the M3Y-type and the DIS 
interactions. As pointed out in Ref. [l3|, difference among the saturating forces is not large at p < po- At relatively 
high density (p > 0.3 fm~^), the M3Y-P3 to P5 interactions have higher E than M3Y-P2 and DIS, reflecting higher 
/C. f (p) of M3Y-P4 and P5 is close to each other even at p w 0.6fm~'^(« 4po)- 

In Fig. [21 contributions of the SE, TE, SO and TO channels in v^^^ -|- v^^^ to E of the symmetric matter are 
shown as a function of fcp. The contribution of the TE and the SO channels in M3Y-P4 is hard to be distinguished 
from that in M3Y-P5. So is the contribution of the TE channel in M3Y-P2. The TE channel takes a minimum at 
fcp ~ 1.5 fm"^, primarily responsible for the saturation at /cp = ^fo ~ 1.3 fm~^. Both the SO and the TO channels 
do not contribute to E significantly at p < po (*-e- ^ ^fo)- While the SO channel becomes attractive in the DIS 
interaction, it is repulsive in the M3Y-type interactions at p > po- The TO channel is repulsive in M3Y-P3 and P5, 
while attractive in M3Y-P4, although the attraction in M3Y-P4 is not so strong as to cause spin polarization in the 
pure neutron matter up to p = 1.3fm~'^(« 8po). Remember that the odd channels in M3Y-P3 are unchanged from 
M3Y-P0. 

Energy per nucleon in the spin-saturated neutron matter (i.e. r\t — —1) is presented in Fig. [31 The result from 
a microscopic calculation in Ref. [l^l is also shown as a reference. The unphysical behavior at high p in the DIS 
result, which comes from the absence of density-dependence in the SE channel, was pointed out in Refs. [l3)[l3- The 
present M3Y-P3 to P5 interactions have relatively strong p-dependence at high p for the neutron matter, if compared 
to M3Y-P2 and DIS. This originates in a(s^)(= 1), and tends to make E{p) closer to the result of Ref. [13. 

IV. PROPERTIES OF DOUBLY MAGIC NUCLEI 

We next turn to doubly magic nuclei, for which the spherical HE approach is expected to be a good approximation. 



6 



TABLE III: Landau-Migdal parameters at the saturation point. 





DIS 


M3Y-P2 


M3Y-P3 


M3Y-P4 


M3Y-P5 


/o 


-0.369 


-0.357 


-0.276 


-0.300 


-0.336 


/l 


-0.909 


-1.044 


-1.027 


-1.005 


-1.112 


/2 


-0.558 


-0.436 


-0.355 


-0.429 


-0.367 


/3 


-0.157 


-0.210 


-0.184 


-0.210 


-0.182 


/o 


0.743 


0.607 


0.578 


0.538 


0.502 


fi 


0.470 


0.635 


0.670 


0.797 


0.692 


fi 


0.342 


0.245 


0.271 


0.286 


0.267 


/a 


0.100 


0.096 


0.104 


0.106 


0.100 


90 


0.466 


0.113 


0.070 


-0.164 


-0.007 


51 


-0.184 


0.273 


0.214 


0.374 


0.299 


32 


0.245 


0.162 


0.160 


0.190 


0.178 


S3 


0.091 


0.078 


0.079 


0.085 


0.081 


So 


0.631 


1.006 


0.933 


1.136 


1.081 


si 


0.610 


0.202 


0.213 


0.109 


0.087 


S2 


-0.038 


0.040 


0.063 


0.016 


0.029 


53 


-0.036 


-0.002 


0.005 


-0.008 


-0.002 



To all the following calculations of finite nuclei, we apply the recently developed algorithm based on the Gaussian 
expansion method (GEM) In this method we employ the s.p. bases of 

f.iMr) - R,e,ir)[Y^'\r)xa]^i^ ; Rutj{r) = M,,, /exp(-z.r2) , (12) 

apart from the isospin index. Here K ^^-^ (f ) expresses the spherical harmonics and Xa the spin wave function. The 
parameter v — v-^ + ivi indicates a complex number corresponding to the range of the Gaussian. Irrespective to 
nuclide, we adopt the following basis parameters [T^ : 

r i/i = (n = 0,l,--- ,5) 
-r = .o6-^", (n^O.1,2) ' (13) 

k f, 2 

with Vi) = (2.40 fm)^^ and h = 1.25 for each {(.,])■ It is notable that, without parameters specific to mass number or 
nuclide, a single set of the GEM bases is applicable to wide range of the nuclear mass table The Hamiltonian 
is H — Hn + Vc — Hem., where Vc and Hem. represent the Coulomb interaction and the center-of-mass (cm.) 
Hamiltonian, while Hn has been given in Eq. ([T]). The exchange term of Vc is treated exactly, in the same manner 
as the nuclear force Vn- Both the one- and the two-body terms of Hem. are subtracted before iteration. 

The calculated binding energies and rms matter radii of several doubly magic nuclei are displayed in Table IIVI 
The results of the new semi-realistic interactions are compared with those of DIS and M3Y-P2 as well as with the 
experimental data. Infiuence of the cm. motion on the matter radii is subtracted in a similar manner to the cm. 
energies The binding energies of these nuclei obtained from DIS (M3Y-P2) are in agreement with the measured 
values within the 3 MeV (5 MeV) accuracy. Though the accuracy is slightly worse, the new interactions also reproduce 
the binding energies moderately well. M3Y-P3 yields underbinding by about 9 - 17 MeV except for ^o^Pb. For M3Y- 
P4 and P5, maximum deviation in the binding energies shown in Table HVl is ^ 7 MeV. Since correlations due to the 
residual interaction could influence, we do not take this deviation seriously at the present stage. The rms matter radii 
of these nuclei calculated from the semi-realistic interactions are comparable to those from the DIS interaction, in 
fair agreement with the data. 

The s.p. levels in ^°*Pb are depicted in Fig. U) The levels obtained from M3Y-P5 are compared with those from 
DIS and the experimental levels. M3Y-P3 and P4 give s.p. levels similar to, though slightly different from, those of 
M3Y-P5. The experimental s.p. energies are obtained from the levels of the neighboring nuclei; 207,209p|-|^ 207rpj ^^^^^ 
^°^Bi. In the HF results, the overall level spacing is relevant to Mq shown in Table HIl In the usual HF calculations 
we have larger level spacing than in the data, and it is not (should not be) remedied until correlations due to the 
residual interaction (or the w-mass) are taken into account |24l | . We thus confirm that the present interactions yield 
as plausible s.p. levels as DIS does. 
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-20- 



-40- 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 

p [fin-2] 

FIG. 1: Energy per nucleon S = E/A in the symmetric nuclear matter for several efltective interactions. Purple dashed, green 
dot-dashed and red solid lines represent the results with the M3Y-P3, P4 and P5 interactions, respectively. Those with the 
M3Y-P2 and DIS interactions are also displayed for comparison by orange and blue solid lines. 

TABLE IV: Binding energies and rms matter radii of several doubly magic nuclei. Experimental data are taken from Refs. [2^ . 









Exp. 


DIS 


M3Y-P2 


M3Y-P3 


M3Y-P4 


M3Y-P5 


16 


-E 


(MeV) 


127.6 


129.5 


127.2 


118.6 


126.3 


126.1 






(fm) 


2.61 


2.61 


2.61 


2.65 


2.60 


2.59 


24q 


-E 


(MeV) 


168.5 


168.6 


165.7 


158.2 


164.0 


166.7 




^/¥) 


(fm) 


3.19 


3.01 


3.06 


3.08 


3.04 


3.03 


«'Ca 


-E 


(MeV) 


342.1 


344.6 


338.8 


325.2 


337.0 


335.1 






(fm) 


3.47 


3.37 


3.38 


3.42 


3.37 


3.37 


4»Ca 


~E 


(MeV) 


416.0 


416.8 


411.9 


401.0 


409.4 


414.1 




^/¥) 


(fm) 


3.57 


3.51 


3.53 


3.56 


3.52 


3.50 


^°Zr 


-E 


(MeV) 


783.9 


785.9 


779.4 


767.9 


775.1 


779.8 






(fm) 


4.32 


4.24 


4.25 


4.27 


4.24 


4.23 


"^Sn 


-E 


(MeV) 


1102.9 


1104.1 


1099.0 


1089.3 


1095.7 


1098.4 






(fm) 




4.77 


4.79 


4.81 


4.77 


4.76 


208pj3 


-E 


(MeV) 


1636.4 


1639.0 


1636.5 


1635.2 


1632.1 


1633.2 






(fm) 


5.49 


5.51 


5.54 


5.55 


5.51 


5.51 



In Ref. it has been shown that the Z — N — 20 shell gaps are narrowed by the tensor force. It is also true in 
the M3Y-type interactions. The s.p. energy difference er(0/7/2) — £T{0d^/2) [t — P, n) in ^"Ca obtained by the DIS 
interaction is in good agreement with the experimental values, both for protons and neutrons, as viewed in Fig. [S] 
While M3Y-P2 and P4 give almost the same size of the shell gaps as DIS does, we have narrower gaps in the HF 
calculations with M3Y-P3 and P5. However, the shell gaps do not collapse by wC^^) of M3Y, in contrast to the 
zero-range tensor force of Ref. We still have 5.2 MeV (7.7 MeV) gap for the proton (neutron) orbits with M3Y-P5. 
These gaps are close to those obtained from the tensor-free Skyrme interaction 'T22' in Ref. It is also worth 
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FIG. 3: Energy per nucleon £. — E/A in the neutron matter for several effective interactions. Circles are the results of Ref. [27|. 
See Fig. [1] for the other conventions. 
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FIG. 4: Single-particle energies for ^''^Pb. Experimental values are extracted from Refs. pjll. |32]. 

commenting that, for M3Y-P5, the octupole correlations significantly influence the ground state of ''°Ca, as will be 
discussed elsewhere. This can make the shell gap look wider, having possibility to account for the observed gap. 

V. PAIRING PROPERTIES 



The M3Y-P2 interaction seems to have reasonable characters in the HF regime, as exemplified in Table HVl However, 
M3Y-P2 has too strong pair correlations, indicating too strong attraction in the SE channel at low densities though 
almost invisible in Fig. [2l This character is inherited from the original M3Y interaction. We have developed the M3Y- 
P3 to P5 parameter-sets by taking the pairing properties into consideration. In this section we shall show characters 
of the new interactions with respect to the pairing. We restrict ourselves to the pairing among like-nuclcons, as usual. 

We implement the spherical HFB calculations for finite nuclei, using the GEM bases of Eqs. (|12I13[) together 
with the £ < 7 truncation. The blocked HFB calculations are applied to the odd-mass nuclei, by assuming that a 
quasiparticle (q.p.) occupies a specified spherical orbital. When several q.p. levels lie closely in energy, we compare 
the total energies by filling each q.p. level and adopt the lowest-energy solution. 



A. Even-odd mass difference in Sn isotopes 

The t^i^^ {n = 1, 2) and tp^^' parameters of M3Y-P3 to P5 are adjusted to the even-odd mass differences of the 
Sn isotopes with 66 < iV < 80. For the mass difference we use the three-point formula A^^^^{N) — E{Z,N) — 
\ [E{Z, N + E{Z, iV - 1)] , with Z = 50 and = odd. The mass differences calculated with M3Y-P4 and P5 
are displayed in Fig. [51 in comparison with the experimental data and with those of DIS. Though not shown to keep 
the figure viewable, M3Y-P3 gives similar A^^gg°(iV) to M3Y-P4. The calculations are not fully convergent for the £ 
truncation. Moreover, the restoration of the particle-number conservation [33] and the non-spherical mean fields [s^l 
could influence A^^gg(A^). Each of them could vary A^^gg{N) by up to a few hundred keV. Not attempting fine 
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FIG. 5: Single-particle energies for Ca. Experimental values are extracted from Refs. |3ll, |32 



tuning of the parameters, we just point out that some of these effects tend to compensate one another in the mass 
differences, and that the new interactions give A^^'^^^{N) to comparable accuracy to the DIS interaction in the same 
model space. 

At ~ 50, 64 and 90, we find that the calculated mass differences depend on the interactions. This is ascribed to 
interaction-dependence of the shell structure. Irregularity and discrepancy at = 63 and 65 should be relevant to 
the A^ = 64 subshell. In M3Y-P5 the subshell effect seems stronger than in the other interactions. At A^ ~ 90 all the 
M3Y-type interactions yield larger mass differences than DIS. This takes place because nl/7/2 and ri2p3/2 well mix 
due to the pairing, in the M3Y-type interactions. At A^ ~ 50 M3Y-P5 yields larger mass difference than the other 
interactions. This is traced back to appreciable excitation from to nOgf/2, which takes place since these two 

orbits are close in energy. Possibly carrying information of the shell structure, data on the masses in A^ ~ 50 and 
A^ ^ 90 will be of interest. 



B. Pairing gap in nuclear matter 



We next view the pairing property in the nuclear matter obtained from the new semi-realistic interactions. In 
phenomenological studies using the Skyrme energy density functionals, it has been argued (ssj whether and how 
much the pair correlations are dominated by the nuclear surface region. Results of the semi-realistic interactions for 
the nuclear matter may provide certain information on this point. However, we find that, in calculating the pairing 
properties, the Yukawa function gives quite slow convergence for the maximum momentum of the s.p. states. In 
practice, even if we cut off the momentum at fc = 50fm~^, which corresponds to £ w 50GeV, the pairing gap is not 
yet fully convergent. It is impractical to include such high energy states in calculations of finite nuclei. As shown in 
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FIG. 6: Even-odd mass differences in the Sn isotopes, A^'^^°{N). The results of DIS, M3Y-P4 and M3Y-P5 are shown by blue 
open circles, green open squares and red diamonds, respectively. Experimental values, presented by black crosses, are taken 
from Ref. [3l[. 



the preceding subsection, we have fixed the SE channel parameters from the HFB calculations of the Sn nuclei using 
the basis-parameters of Eq. (|13p . High momentum components are automatically excluded in the basis set. It will be 
natural to introduce a certain cut-off in arguing the pairing in the nuclear matter, and the cut-off should desirably be 
consistent with the basis set of Eq. (|13p . 

The basis set is composed of radial Gaussians, whose Fourier transforms are again Gaussians in the momentum 
space. We here consider a cut-off factor for the s.p. momentum space of 



g{k) = 0{kc - k) + 0{k - fc^) exp 



(14) 



The measure in the k integration is multiplied by g{k). Among the bases of Eq. (jl3[) . the highest k component is given 
by the i' — 1/0(1 ± ^i) basis. Since the Fourier transform of this basis is proportional to exp[ — fc^(l ± ■|j)/4t'o(l + ^)] , 

k ko = 2^vq{\ + ^){~ 1.55 fm^^) components are well included in the set. To be consistent with the basis set, it 
will be reasonable to take kc ~ fco for the nuclear matter calculation. We here consider three cases, [kc, kj) = (fco, ^o)j 
(2fco, fco) and (4fm~^, 0). The last choice of the fc^ — > limit indicates a sharp cut-off, and k^ = 4fm^^ approximately 
corresponds to the maximum q.p. energy in the HFB calculations of the Sn nuclei. 

In Fig. [71 the pairing gap at the Fermi energy, which is obtained from the Bardeen-Cooper-Schrieffer (BCS) 
calculation in the symmetric nuclear matter using the method of Ref. [s^ , is plotted as a function offcp ^ {iir'^ p/2f/^ . 
The pairing among like-nucleons arises from the SE channel in the effective interaction. The result of the M3Y-P5 
interaction is compared with that of the Gogny DIS interaction. The cut-off is not needed for the Gogny interaction, 
and the present cut-off does not influence the DIS gap. The pairing gaps of M3Y-P3 and P4 are similar to that of 
M3Y-P5. 

The new semi-realistic interactions are not drastically different from DIS, in respect to the pairing properties in the 
nuclear matter. The gap has a peak at k-p « 0.8 fm^^ « O.Gfcpoj namely at p sa 0.2/3o, for all cases. However, the peak 
height and the behavior at p > 0.2po are different between the M3Y-type interactions and DIS. With rapid decrease 
at jO > 0.2/50, the M3Y-type interactions have more surface-dominant pairing than DIS. In DIS, the pair correlation 
may have sizable contribution from the bulk, even though it is stronger at the nuclear surface. Note that the cut-off 
parameter does not influence the nuclear matter pairing qualitatively, as long as it is more or less harmonious with 
the basis-set adopted in the calculations of finite nuclei. 
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FIG. 7: Pairing gap at the Fermi level ep in the symmetric nuclear matter. Red dot-dashed, solid and dotted curves are 
obtained from M3Y-P5 with the momentum cut-ofT {kc,kd) = (fco,fco), (2fco,fco) and (4fm^^,0) in Eq. (|14|l . respectively. Blue 
solid curve displays gap of the DIS interaction. 

VI. NEUTRON DRIP LINE 

Prediction of the neutron drip line depends on effective interactions to a certain degree. In this section we compare 
location of the neutron drip line predicted by the spherical HFB calculations with the present semi- realistic interactions 
and with the Gogny DIS interaction for the O, Ca and Ni isotopes. Although complete description of the drip line 
may require fine tuning of the parameters as well as taking account of correlation effects, it will be interesting to see 
what is relevant to location of the drip line. 

A. Z = 8 nuclei 

We present the two-neutron separation energies S2n for the O isotopes in Fig.[8l The calculated values are compared 
with the experimental data. We do not show S2n if the neutron chemical potential is positive. Though not displayed in 
Fig. [HI obtained from M3Y-P3 is close either to that from DIS or M3Y-P4. Whereas 25-28q j^^ve experimentally 
been established to be unbound |37|], most MF calculations so far have failed to reproduce this nature. It would be 
noteworthy that one of the present semi-realistic interactions, M3Y-P5, correctly describes the location of the neutron 
drip line for oxygen within the spherical HFB approximation; ^^O is the heaviest bound oxygen isotope. The ^^O 
nucleus has higher energy than ^'^O, and in 26-28q ^j-^g chemical potential becomes positive. In contrast, ^^O is bound 
in the HFB calculations with M3Y-P3 and P4, as in the calculation with the DIS interaction. 

The reason why ^^O is not bound with M3Y-P5 can be traced back to £71(0^3/2), the s.p. energy of nO(i3/2. In 
Fig. [51 the neutron s.p. energies in the HF calculations are depicted. We show the energies obtained from the bases of 
Eq. (jl3p even when the s.p. energy is positive, not treating the boundary condition carefully. We view that M3Y-P5 
gives higher £71(0^3/2) than the other interactions, which originates from the slightly stronger u^'"^) as well as from 
the relatively small Mq . Note that we have fixed the enhancement factor for v^^^'> in M3Y-P5 so as to reproduce the 
s.p. spectrum around ^°^Pb, not adjusting accurately to, e.g., the is splitting around ^^O. It is commented that 
has small but attractive contribution to £71(0^3/2), and therefore is irrelevant to the higher £„(0(i3/2) in the M3Y-P5 
result. 
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FIG. 8: 5*2,1 of the O isotopes (A'' — even). The calculations are performed in the spherical HFB approximation. Experimental 
data are taken from Ref. (Sid. See Fig. [6] for conventions. 
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FIG. 9: HF single-particle energies in the O isotopes. Blue, green and red lines represent the results with the DIS, M3Y-P4 and 
M3Y-P5 interactions, respectively. For each interaction, dot-dashed line is for en(0d5/2)) solid line for £n(lsi/2) and dashed 
line for e„ (00(3/2). 
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TABLE V: Neutron numbers of the heaviest bound Ca and Ni nuclei predicted by the spherical HFB calculations with several 
interactions. 



Isotope 


DIS 


M3Y-P3 


M3Y-P4 


M3Y-P5 


Ca 


44 


50 


48 


50 


Ni 


58 


64 


62 


60 




FIG. 10: Difference between the HF and HFB energies for the Ca isotopes (A'' — even), obtained from DIS, M3Y-P4 and P5. 
See Fig. [6] for conventions. Dotted lines are drawn to guide eyes. 



B. Z = 20 and 28 nuclei 



Location of the neutron drip line for the Ca and Ni nuclei could be investigated by the currently constructed or 
designed experimental facilities [38| . We tabulate location of the neutron drip line predicted by the spherical HFB 
calculations with the M3Y-type and the DIS interactions, in Table fVl 

If we use the DIS interaction, the heaviest bound Ca nucleus is ^''Ca, because the neutron chemical potential is 
positive in > 45. The M3Y-P3 and P5 interactions predict that ^"Ca is bound, while ^^Ca is the heaviest bound 
Ca isotope in the M3Y-P4 result. We depict difference between the HF and the HFB energies for the Ca nuclei in 
Fig. [ini which represents the pair correlation. Though not shown, M3Y-P3 gives similar results to M3Y-P5. While 
the pairing effects are in good agreement among all of these interactions in A^ < 32, the M3Y-type interactions give 
stronger pairing than DIS in A^ > 34. 

Both in the predicted position of the neutron drip line and in the pair correlation in A^ > 34, the s.p. energy of 
nOf/g/2 plays an important role. The neutron drip line can extend up to ™Ca if £71(039/2) is sufficiently low. At ^"Ca 
we have 6^(099/2) = +0.73, -0.63, +0.23 and -0.65 MeV in the HF calculations with DIS, M3Y-P3, P4 and P5, 
respectively, well correlated to the location of the drip line. The lower en{0g9/2) leads to the smaller shell gap at 
A^ = 40, Sn{0g9/2) — £11(0/5/2), which makes the pair excitation across A^ = 40 easier. The gap is 3.8, 2.0, 3.7 and 
1.7 MeV in DIS, M3Y-P3, P4 and P5. It is noted that this quenching of the shell gap in M3Y-P3 and P5 comes from 
the relatively strong ti^^^) to some degree, but not from u^^^), because u^^^) hardly contributes to the s.p. energies 
in an £s-closed shell. 

The vanishing difference between the HF and HFB energies is often connected to the shell (or sub-shell) closure. 
FigurefTPlindicates that, while ^°Ca is stiff against the pair excitation with DIS and M3Y-P4, significant pair excitation 
occurs by M3Y-P5 (and by M3Y-P3), because of the small shell gap at A^ = 40. On the contrary, A^ = 50 is stiff against 
the pair excitation in the present M3Y-P3 and P5 results, in which '^°Ca is bound. We have £^(0^9/2) ~ —2 MeV 
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FIG. 11: Difference between the HF and HFB energies for the Ni isotopes (A'^ = even). See Fig. |6]for conventions. 



at Ca, and the gap between nOgg/2 and the continuum seems large enough for iV = 50 to keep the magic nature 
against the pairing, if we use the present semi-reahstic pairing interaction. Figure [10] also suggests shell closure at 
N = 32. We shaU return to this point in Sec lVlIl 

In Fig. [TTl difference between the HF and HFB energies is shown for the Ni isotopes. For all the interactions, the 
energy difference becomes vanishingly small at N = 20, 28, 40, 50 and 58. We view significant interaction-dependence 
in 28 < iV < 40. In particular, the pair correlation is suppressed at iV = 32 with DIS, while no such effect is found 
with M3Y-P5. It is noted that, with M3Y-P5, ^'^Ni seems almost doubly magic as is consistent with experiments [s^, 
although ^"Ca is not, as has been seen in Fig.fTUl For the neutron-rich Ni region, energy sequence of the s.p. orbitals 
above iV 50 is 2si/2, 1^3/2 and 0177/2, from the lower orbit to the higher. The hindrance of the pair excitation 

at ®^Ni suggests magic or submagic nature of = 58 due to the gap between n2si/2 and n\d^/2- Unlike the M3Y- 
type interactions, by DIS the pair excitation is hindered also at A'^ = 56. The predicted neutron drip line appreciably 
depends on en{^d^/2) and £n{^g7/2)- The DIS interaction yields higher £71(1^3/2) than the M3Y-type interactions, 
which causes the drip line at iV = 58. The low Sni^gi /2) in M3Y-P3 induces pair excitation to 037/2, leading to the 
binding up to ^^Ni. 

In the highly neutron-rich region, the diffuseness of the nuclear surface becomes larger than in the /^-stable region. 
Then the pair correlation could be relatively strong if it has the surface-dominant nature. However, it is not easy 
to argue precisely the extent of surface dominance in the pairing from location of the drip line, since it is obscured 
by influence of the shell structure. In practice, when we use M3Y-P5 for the HF Hamiltonian and DIS for the pair 
potential, predicted location of the drip line for Ca is the same as the result of the pure M3Y-P5 prediction. For 
the Ni case, although *^Ni becomes unbound when DIS is used for the pair potential, the chemical potential is only 
— O.OTMeV in the pure M3Y-P5 result. It will be fair to say that the difference in the pairing channel between the 
M3Y-type interactions and DIS is not quite significant to location of the neutron drip line. Dependence of the rms 
matter radii on the pairing interaction is not apparent either, as long as we work with the present M3Y-type or the 
DIS interactions, whereas the radii near the drip line are sensitive to the separation energies. 



VII. NUCLEUS-DEPENDENCE OF SINGLE PARTICLE ENERGIES 



It has been pointed out that the shell structure, particularly its nucleus-dependence (sometimes called shell evolu- 
tion), may be connected to characters of effective interactions [1, [l^, In recent studies role of the tensor force in 
the shell structure has attracted great interest [1, Hi ■ It is known that observed energies of one-particle states on top 
of a certain core are appreciably disturbed by correlations beyond the MF regime. However, it is presumable that 
those correlations do not vary in a certain region of nuclei. In the Sb and the A'^ = 83 nuclei, s.p. energies of a few 
orbitals are extracted from several fragmented states [41| , by averaging their energies weighted by the spectroscopic 
factors. As a result, the averaged s.p. energies are shifted from the lowest states with specific spin-parity nearly by a 
constant, from nucleus to nucleus. While data on the averaged energies are not available in many cases, we proceed 
to investigate nucleus-dependence of the shell structure by using the measured energies of the lowest states. 
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FIG. 12: Asn ~ en(0d3/2) — £n(lsi/2) for the A'^ = 16 isotones. Blue, green and red lines correspond to the results with the 
DIS, M3Y-P4 and P5 interactions, respectively. Thin red solid and dashed lines represent relative contributions of v^™'' and 
"np'pp in the M3Y-P5 resuh. 



In this section we shall investigate nucleus-dependence of s.p. energies, using the spherical HF or HFB calculations. 
Results of several interactions are compared. In the HP calculations for open-shell nuclei, the HF Haniiltonian is 
obtained by folding the interaction by the occupation numbers on each spherical orbital up to the Fermi level. 



A. Neutron orbits and shell gap in A'^ = 16 nuclei 

Nucleus-dependence of s.p. energies could be relevant to the new magic numbers in unstable nuclei [3, IS]- We 
investigated the s.p. energies in the TV = 16 and 32 nuclei in Refs. [13, El, and disclosed role of Wq^^ep for the = 16 
isotones, using the M3Y-P2 interaction. We shall reinvestigate nucleus-dependence of the s.p. energies in these nuclei, 
drawing attention also to the tensor force. 

^'-dependence of the neutron s.p. energy en{0d^/2) relative to en(lsi/2) is appreciably affected by effective in- 
teractions [13, US]- Figure [T2l depicts Z\£„ = en(0(i3/2) — en(lsi/2) for varying Z obtained from the HF calcu- 
lations in the TV = 16 isotones. Though not shown, Z\£„ of M3Y-P3 resembles that of M3Y-P5. To clarify 
role of wC^^) and Wqpep, we also plot their contributions to Z\e„ in the M3Y-P5 result, which are calculated as 
T,f{Nj') {2J +l){jj'J\v\jj'J)/{2j + l){2j' + 1) and shifted by the values at iV = 14. 

The present semi-realistic interactions yield increasing Asn as Z goes from 14 to 8, in contrast to the DIS interaction. 
We view in Fig.[T2]that w[^^ep B2| and v'^™'^ produce this feature cooperatively. This Z-dependence of the s.p. energies 
could be relevant to the new magic number TV = 16 in the neutron-rich region (43] . It has been confirmed that several 
popular Skyrme interactions show similar behavior to DIS [42| . 
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FIG. 13: Z\e„(0/5/2) = £71(0/5/2) — £n(lp3/2) (solid lines) and Z\£„(lpi/2) = £n(lpi/2) - £n(lp3/2) (dot-dashed lines) for the 
A'' = 32 isotones. See Fig. [T2]for conventions of colors. Thin red solid and dashed lines represent relative contributions of v^™^ 
and vl^p^p to Z\£„(0/5/2) in the M3Y-P5 result. 



B. Neutron orbits and shell gap in A'^ = 32 nuclei 

In Fig.[Tni Z-dependence of the neutron s.p. energies relative to en(lp3/2), ^^nij) = £n(j) —£71(1^3/2), is shown for 
the = 32 nuclei, by taking j — O/5/2 and lpi/2- As in the preceding subsection, contributions of yC^'^) and Wq^^ep 
to Z\£„ (0/5/2) in the M3Y-P5 result are also presented, after shifting by the values at Z = 28. With the M3Y-P5 
interaction we obtain strong Z-dependence in Z\e„(0/5/2), which could be relevant to the magicity of = 32 in the 

neutron-rich region Once again this Z-dependence originates in Vqpep [li| and v^'^^K 

This behavior of ^£71(0/5/2) is reflected in the pair correlations shown in Figs. [TOl and [TT] The pairing effects are 
small at iV = 32 for all the interactions in the calcium case, because Z\£„(0/5/2) as well as Z\e„(lp3/2) are greater 
than 2 MeV. Recall that this s.p. energy difference competes with the pairing gap, whose typical value is estimated 
to be A w 12yl~^/^ « 1.7 MeV. On the other hand, in the nickel case the narrow Z\e„(0/5/2) leads to substantial pair 
excitation at A^ = 32 for the M3Y-type interactions, while such excitation is kept suppressed in the DIS result. 

In the neutron-rich region of Z ~ 20, there was a prediction that A^ = 34 should be a magic number, based on 
a shell model calculation [4^. The present MF calculations with the semi- realistic interactions do not support this 
prediction. While the TV = 32 shell gap is 2.7 MeV for ^^Ca in the HF calculation with M3Y-P5, the iV = 34 gap is 
only 1.2 MeV for ^"^Ca. The pair excitation across A^ = 34 is sizable, as has been viewed in Fig. [TOl 



C. Proton orbits in Z = 50 nuclei 



In recent studies, nucleus -dep endence of s.p. energies in the Sn isotopes and in the = 82 isotones has been 
disclosed from experiments [4l|. It has been pointed out that the tensor force seems to play a crucial role in the 
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FIG. 14: 5Aep{j) in the Sn isotopes (A'^ — even), for j — Og-j-/2 (dot-dashed hnes) and Oftii/2 (sohd Unes). Blue, green and red 
Unes represent the results of DIS, M3Y-P4 and M3Y-P5, as before. Pluses (j — Og-j-/2) and crosses (j — O/111/2) are experimental 
values taken from the lowest states of the Sb nuclei [3^ . Thin red lines are contributions of v'^™^ in the s.p. levels of M3Y-P5. 



TV-dependence of the proton s.p. energies in the Sn isotopes [18|. We now have the M3Y-type interactions with quite 
realistic tensor force (M3Y-P3 and P5) and without tensor force (M3Y-P4), both of which reproduce the properties 
of doubly magic nuclei as well as the pairing properties to reasonable accuracy. We apply the HFB calculations with 
these new interactions to investigating the nucleus-dependence of the s.p. energies in the Sn nuclei in this subsection, 
and in the = 50 and iV = 82 nuclei in the subsequent subsection. 

In the previous studies [iMIillj the relative proton s.p. energies ep(0/iii/2) —£piOg7/2) were the point of discussion. 
We here consider the energies of these two orbits relative to 1^5/2, ^£p{j) = £pU) — ep(l'^5/2) with j — O57/2 and 
0/iii/2- Taking Aep{j) at TV = 64 to be a reference and denoting it by Aep{j), we plot SAep{j) = Aep{j) — Aep{j) 
in Fig. [131 The values of the M3Y-type interactions (M3Y-P4 and P5) are presented together with those of DIS and 
the experimental data. For the data we use the energies of the lowest states in the Sb nuclei. 

The semi-realistic M3Y-P5 interaction reproduces 5Aep{j) remarkably well. In the 64 < N < 82 region, the 
A'^-dependence of Aep{j) takes place due primarily to the occupation of nO/iii/2, to which contribution of 
significant. Though not shown, 5Aep{j) with M3Y-P3 is close to the M3Y-P5 result. With DIS, whereas 5Aep{0gi/2) 
is in qualitative agreement with the data, SAep{Qhiif2) is not, because the tensor force is absent. The same holds 
for M3Y-P4. Although M3Y-P4 reproduces the tendency of the observed iV-dependence of ep(0/iii/2) — £p(057/2), 
it gives wrong behavior for SAep{0hii/2)', i-e., iV-dependence of ep(0/iii/2) relative to ep(ld5/2)- It is also remarked 
that M3Y-P5 gives quite different behavior of 5Aep{j) from DIS and M3Y-P4 in iV < 64, and that the currently 

(c) 

available data favor the result of M3Y-P5. We have confirmed that Wopep does not have important effects on the 
TV-dependence of Aep{j). 

It is emphasized that the M3Y-P5 interaction can reproduce the variation of the s.p. energy difference in the 
Sn isotopes without destroying the shell structure of the doubly magic nuclei shown in Sec. IIVI This could be an 
advantage of the present realistic tensor force. In Ref. ^4||, no such parameters were found within the Skyrme density 
functional including the zero-range tensor force. 
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FIG. 15: SAsnij) ™ the A'^ = 50 isotones {Z — even), for j = OgY/2 and 0/iii/2- Conventions are the same as in Fig. 1141 except 
that j represents the neutron orbits and Z = 40 is taken to be a reference. Experimental values are taken from the lowest 
states of the A'^ = 51 nuclei [3^ . 



D. Neutron orbits in iV = 50 and N = 82 nuclei 

In the preceding subsection, we have seen that the tensor force affects Sp{0gi^2) ^^^d £p(0/iii/2) via the occupation 
of nO/iii/2- This is accounted for by the attractive (repulsive) nature of the tensor force between a neutron occupying 
a j> — i + 1/2 orbit and a proton occupying = £' — 1/2 {j'^ = i' + 1/2) [l3|. The same mechanism is expected 
for £71(057/2) and e„(0/iii/2) in the = 50 nuclei, as pOgg/2 is occupied. We define Z\£„(j) — e„(j) — e„(ld5/2), 
and take its value at Z = 40 to be Z\e°(j'). In Fig. [15] JZ\e„(j') — Z\£„(j) — Z\e°(j) is displayed for j = 057/2 and 
O/111/2. Although £71(057/2) varies almost in parallel to £n(0/iii/2) in Z > 40 for DIS, notable Z-dependence arises 
for M3Y-P5, in qualitative agreement with the observed s.p. levels. The tensor force has significant contribution to 
this behavior as in the Z — 50 case. More realistic than DIS in the central channels but not having the tensor force, 
M3Y-P4 yields (5Z\£„(j) in-between which is qualitatively good but quantitatively insufficient. 

For the = 82 nuclei, we consider Z\£„(j) — £„(j) — £71(1/7/2), from which SAsnU) — Aen{j) — Zi£°(j) is obtained 
by assuming the value at Z — 64 as Z\£°(j). Figure [TBI shows SAsnij) fo'" J — O/ig/2 and 0ii3/2. The M3Y-P5 
interaction well describes (5Z\£„(j) in Z > 64, which is affected mainly by the occupation of pO/iii/2. However, we 
cannot fully reproduce the tendency in Z < 64, which was argued in Ref. [4l|, although the M3Y-P5 results of dAsnij) 
are substantially better than those of DIS and M3Y-P4. Further investigation will be necessary. 

VIII. SUMMARY AND OUTLOOK 

We have developed semi-realistic effective interactions to describe low energy phenomena of nuclei. Starting from the 
M3Y interaction, we add a density-dependent contact force and modify several strength parameters in a phenomeno- 
logical manner, whereas maintaining the OPEP part in the central force. We have obtained three new parameter-sets; 
two of them (M3Y-P3 and P5) keep the tensor force of the M3Y-Paris interaction, and the other (M3Y-P4) has no 
tensor force. Basic characters of the interactions are checked by the Hartree-Fock calculations for the infinite nuclear 
matter, and for the doubly magic nuclei. The singlet-even channels of the interactions, which are relevant to the pair- 
ing properties, are fixed from the even-odd mass differences in the Sn isotopes, by using the Hartree-Fock-Bogolyubov 
calculations. 

We further implement the Hartree-Fock and the Hartree-Fock-Bogolyubov calculations for spherical nuclei, applying 
the new interactions. Predicted shell structure depends on the effective interactions to certain extent. This may 
significantly affect location of the drip lines. The new semi-realistic interaction M3Y-P5 correctly describes the 
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FIG. 16: SAen{j) in the N = 82 isotones {Z = even), for j = 0/19/2 (dot-dashed lines) and 0ii3/2 (solid lines). See Fig. 1141 for 
conventions of colors. Pluses {j = 0/19/2) and crosses (j = 0113/2) are experimental values taken from the lowest states of the 
iV = 83 nuclei fsl. 



experimental consequence that the heaviest bound oxygen is ^ O. We have argued location of the neutron drip line 
for the Ca and the Ni nuclei, and its relevance to the shell structure. Variation of the single-particle (s.p.) energies, 
particularly contribution of the tensor force to it, is a current topic. It is suggested that the tensor force as well as 
the OPEP part of the central force play a significant role in the magic numbers A'^ = 16 and 32 in the neutron-rich 
region. It has been shown that the new semi-realistic interactions including the tensor force, M3Y-P5 in particular, 
describe the variation of the s.p. levels fairly well in Z = 50, = 50 and iV = 82 nuclei. It is remarked that this 
interaction can reproduce the variation of the s.p. energy difference in the Sn isotopes without destroying the shell 
structure of the doubly magic nuclei. 

It will be of interest to apply the semi-realistic interactions to deformed nuclei, and to excited states via the 
random-phase approximation (RPA). Both projects are in progress (for the latter, see Ref. [1^). 
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